Super-resolution ghost imaging via compressive sampling reconstruction 
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For ghost imaging, pursuing high resolution images and short acquisition times required for re- 
constructing images are always two main goals. We report an image reconstruction algorithm called 
compressive sampling (CS) reconstruction to recover ghost images. By CS reconstruction, ghost 
imaging with both super-resolution and a good signal-to-noise ratio can be obtained via short ac- 
quisition times. Both effect influencing and approaches further improving the resolution of ghost 
images via CS reconstruction, relationship between ghost imaging and CS theory are also discussed. 
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In recent ten years, ghost imaging (GI) has attracted 
lots of attentions in the field of quantum optics [I] . The 
image of an unknown object can be nonlocally recon- 
structed by the intensity correlation measurements be- 
tween two light fields. Both entangled source and ther- 
mal light can be used to realizeghost im aging [1 H H 
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researches of ghost imaging have demonstrated that all 
the information of an object can be obtained by "global 
random" measurements and lots of problems which are 
hard to be solved by conventional imaging approaches, 
such as x-ray diffraction imaging, projection imaging in 
Fraunhofer region and imaging in scatte ring media, can 
be settled by ghost imaging H El El E2j- However, 
the reconstruction algorithm of intensity correlation mea- 
surements (also called GI reconstruction) always faces 
with two main drawbacks. One is that the best resolu- 
tion of recovered images is determined by the size of the 
speckle placed on the object plane based on the previ- 
ous experimental results [H Ell • The other is that long 
acquisition times should be required for reconstructing 
images with a good signal-to-noise ratio (SNR) [l[ . 

Recently, compressive sampling (CS) theory has 
proved mathematically that "global random" measure- 
ments have much higher image extraction efficiency than 
"point-by-point" scanning measurements [H [H [H [22|, 
HI HIl- Correspondingly, an advanced reconstruction 
algorithm called CS reconstruction has asserted experi- 
mentally that one can recover certain signals and images 
from far fewer samples or measurements than traditional 
methods use [H HI [H [H Hi| . CS reconstruction also 
offers great potential for better resolution over classical 
imaging [25| . However, in these experiments, the pseudo- 
thermal incoherent measurement matrix A is obtained by 
hardware methods and the test detector should collect 
all the intensity from the object. In practical imaging 
applications, it is impossible to collect all the intensity 



from the object in most of the imaging schematics, but 
the complete image of an object can still be realized by 
ghost ima ging even if the test detector is a pointlike de- 
tector [H13 Eg]- In this letter, super-resolution ghost 
imaging via CS reconstruction is investigated when the 
test detector just collects partial intensity form the object 
and relationship between ghost imaging and CS theory 
is discussed. 

Fig. 1(a) represents standard schematic of conven- 
tional imaging. The experimental schematic for ghost 
imaging with thermal light is shown in Fig. 1(b). The 
light source S, which is obtained by passing a laser beam 
through a slowly rotating ground glass disk fl2| , first 
propagates through a beam splitter, then is divided into 
a test and a reference path. In the test path, the partial 
intensity which is transmitted through the object is col- 
lected by the test detector D t . In the reference path, a 
copy of the speckle field which impinges on the object is 
recorded with a CCD camera D r . 

By optical coherence theory the intensity distri- 
bution obtained by first-order correlation can be repre- 
sented as: 

I(x) — I dxi / dx 2 G i - 1 ' 1 \xi,x 2 )h* t {x,xi)h t {x,x 2 ). (1) 



where G^ 1,1 ' {x\ ,x 2 ) is the first-order correlation function 
on the source plane. h t {x,xi), h* t {x,x\) are the impulse 
function of optical system and phase conjugate of the 
impulse function, respectively. 

Based on ghost imaging via GI reconstruction [1 [||, 
we can obtain the correlation function between the de- 
tectors: 
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xh*(x 1 ,x r )ht(x 2 ,x t ) | . (2) 

where ht(xt,x 2 ) is the impulse function in the test path 
whereas h*(x r , x±) denotes phase conjugate of the im- 
pulse function in the reference path. 

Suppose the light source is fully spatially incoherent, 
then 
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G^(x 1 ,x 2 )=I S(x 1 -x 2 ). 



(3) 
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FIG. 1: Schematics of conventional imaging and ghost imag- 
ing with thermal light, (a). Conventional imaging; and (b). 
Ghost imaging; (1) and (2) are the schemes of the test path 
for the experiment and the discussion of further improving 
the resolution of ghost images, respectively. 



where Iq is a constant, and S(x) is Dirac delta function. 

For the schematic shown in Fig. l(b-l), under the 
paraxial approximation, the impulse response function 
of the reference system is 



h r (x r , %i) oc exp{ — (x r — xi) 2 }. 
Xz 



(4) 



When -i- + = j, then the impulse response function 
for the test path is 

ht(xt,X2) oc / dx' exp{^-(x' - x 2 ) 2 }t{x') 
J Az 

xexp{-^' 2 }sin C fe + £)]}. ( 5 ) 
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where t(x), L are the transmission function of the object 
and effective aperture of the imaging lens /, respectively. 
Furthermore, sin c(x) = SS^l. Substituting Eqs. (3) 
and (5) into Eq. (1), the intensity distribution on the 
CCD camera in Fig. 1(a) is 



I(x) cx f dx' \t(x')f sinc 2 [v(— + — )]. 
J A z 2 zi 



(6) 



Similarly, substituting Eqs. (3)- (5) into Eq. (2), the 
correlation function for ghost imaging can be represented 
as 



AG^ 2) (x r ,x t ) cx 



dx' sincf— (x r — x')]t(x') 
\z 
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where D is the transverse size of the source. If the small- 
est length scale of the object is larger than the size of 
the diffraction limit cased by the lens' effective aperture 
L (namely the lens / can nearly collect all the intensity 
from the object), thus sinc[y (f t + f^)] ~ ^(fj + ff )> tnen 

AG (2 ' 2) (s r )= / dx t AG^ 2 - 2) (x ri x t ) 



(8) 



dx' |i(x')| 2 sinc 2 [— (x r — x')). 
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From Eqs. (6) and (8), the best resolution of images 
are determined by the size of the speckle placed on the 
object plane (namely Aa; ~ ^j) for ghost imaging via GI 
reconstruction whereas for conventional imaging depend 
on the effective aperture of the imaging lens / (namely 

However, for ghost imaging, except for GI reconstruc- 
tion to recover the object's transmission function, an- 
other novel method called CS reconstruction algorithm 
perfectly accords with the physical principle of ghost 
imagin g an d can be used to recover the object. By CS 
theory [H HI EH 0, HE HI , when the image 
to be recovered with n x n matrix can be viewed as a 

2 

vector in R n , the measurements process can be written 
compactly in matrix notation as 



y = Ax. 



(9) 



where x G R n is the "true" image. The m x R n matrix 
A is constructed by stacking the m measurement basis 
functions — ► each of which is also a vector in R n and 
any two of which are mutually incoherent, y is the m vec- 
tor containing the observations. Given y, we reconstruct 
the image by solving the following convex optimization 
program [20 |: 



min ||jc||. subject to y = Ax. 



(10) 



where ||V|| f = J2i \ v i\ * s the l\ norm of V . 

In the schematic of ghost imaging shown in Fig. 1(b), 
the m x R n matrix A can be obtained by the reference 
path and each of the measurement basis functions is con- 
structed by one of realizations recorded with the CCD 
camera D r . Because each of the speckle fields registered 
by the CCD camera D r is random and independent of 
others, thus the property of incoherence sampling among 
the m vectors of matrix A is obviously satisfied. Simi- 
larly, the m vector y is corresponded to the m measure- 
ments registered by the bucket detector D t . If the speckle 
field which impinges on the object is described by I r (x, y) 
and B r denotes the total intensity recorded by the detec- 
tor D t , then the image of the object can be reconstructed 
by solving the following convex optimization program by 
Eq. (10): 

tcs — t which minimizes : ||t (a;,y)||j ; subject to 

J dx J dyl r (x,y)t'(x,y) = B r ,V r = 1 • • -m.(ll) 
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FIG. 2: Simulated conventional imaging, ghost imaging via 
GI and CS reconstructions of a double-slit with z=200mm, 
/=250mm, zi=Z2=500mm and L=6mm when the test path 
of ghost imaging is proposed as Fig. l(b-l). (a). The ob- 
ject; (b). Conventional imaging; (c). Ghost imaging via GI 
reconstruction with 2000 realizations; (d) and (e) are ghost 
imaging via CS reconstruction with 32 realizations when the 
pixel resolution of the CCD camera D r is 3/im and 18/im, 
respectively. 

where tcs is the object's transmission function recovered 
by CS reconstruction algorithm. 

For CS reconstruction, the resolution of reconstructed 
object depends on the row vector of matrix A (namely 
R n ), while the image of an object is obtained by scan- 
ning the position of the photons on the CCD camera D r 
for ghost imaging. Thus the resolution of ghost imaging 
recovered via CS reconstruction is closely related to the 
resolution of the CCD camera D r to record I r (x, y). 

In the experiment, the wavelength of the source was 
A=650nm and the transverse size of the source was 
D=2.0mm. Fig. 2 presents the numerical simulation re- 
sults for conventional imaging, ghost imaging via GI and 
CS reconstructions when the test path of ghost imaging is 
proposed as Fig. l(b-l). The object is a double-slit with 
slit width a=30/im, slit height /i=120^m and center-to- 
center separation d=60/xm. For the CS reconstruction, 
we have utilized the ^i-Magic algorithm 0, [26j. From 
Eqs. (6) and (8), the resolution of both conventional 
imaging and ghost imaging recovered by GI reconstruc- 
tion are very low, which are also demonstrated in Fig. 
2(b-c). However, a high-resolution ghost image can be 
obtained by CS reconstruction (Fig. 2(d)). Correspond- 
ingly, from Fig. 2(d-e), the resolution of recovered images 
via CS reconstruction will reduce as the resolution of the 
CCD cameras D r is decreased. 

For ghost imaging, the lens / shown in Fig. l(b-l) is 
just used to collect the information from the object and 
the lens' effective aperture L limits the information of 
high frequency. In practical applications, we can use a 
device with large collecting area to collect the informa- 
tion from the object instead of the lens /, which is similar 
to the case in Fig. l(b-2). In Fig. 1(b), by replacing the 
test path (1) using the scheme (2), ghost imaging via GI 
and CS reconstructions of the same double-slit in differ- 
ent collecting areas for the test detector D t are shown in 
Fig. 3. From Fig. 3(a-b), enlarging the collecting areas 
of the test detector can improve the resolution of ghost 
imaging when a < 1 - 2 ^ Zl (namely the collecting areas 
of the test detector is smaller than 13. 2mmx 13.2mm). 
However, similar to the results described in Eq. (8), if 
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FIG. 3: Numerical simulation results of ghost imaging via GI 
and CS reconstructions with zi=500mm when Fig. l(b-2) 
is used as the test path of ghost imaging, (a). The inten- 
sity distribution on the test detector plane with 2000 real- 
izations; (b). GI reconstruction with 2000 realizations; (c) 
and (d), respectively, are CS reconstruction with 32 real- 
izations when the pixel resolution of the CCD camera -Or- 
is 3/xm and 18/xm. The collecting areas of the test de- 
tector Dt with the sizes L\ x L\ shown in (1), (2), (3) 
and (4) are 6mm x 6mm, 10mm x 10mm, 15mm x 15mm and 
30mm x 30mm, respectively. 



a > 1 - 2 ^ Zl , then the best resolution of ghost imaging 
via GI reconstruction will be determined only by the size 
of the speckle placed on the object plane, which is also 
been demonstrated in Fig. 3(b). Correspondingly, more 
information from the object is collected by the test de- 
tector because of the increase of the collecting areas, thus 
the quality of ghost imaging via CS reconstruction will 
be enhanced (Fig. 3 (c-d)). Fig. 3(c-d) also shows that 
the resolution of recovered images via CS reconstruction 
will reduce as the decrease of the resolution of the CCD 
cameras D r . 

For conventional imaging shown in Fig. 1(a), the im- 
age registered by the CCD camera is obtained by long 
exposure. So the intensity distribution of the light field 
on the object plane is uniform and CS reconstruction 
can not be applied because of coherence sampling. Dif- 
ferent from conventional imaging, ghost imaging is based 
on the principle of short exposure and the acquisitions 
recorded with the CCD cameras D r are incoherence sam- 
pling. Thus, the process of ghost imaging shown in Fig. 
1(b) describes vividly the standard CS theory but we also 
demonstrated that ghost imaging can be recovered by CS 
reconstruction even if the intensity from the object are 
partially collected. 

Furthermore, for conventional imaging, enlarging the 
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effective aperture of imaging lens /, the resolution of the 
image can be improved. Also, the resolution of ghost im- 
ages via GI reconstruction can be enhanced by decreas- 
ing the transverse coherence width placed on the object 
plane (such as a source with large transverse size, a small 
distance z). However, the resolution of recovered images 
via CS reconstruction is related to the resolution of the 
CCD camera D r (Fig. 2(d-e), Fig. 3(c-d)). So for ghost 
imaging via CS reconstruction, the resolution of ghost 
images can be further improved by increasing the res- 
olution of the CCD camera D r or computational ghost 
imaging approaches. From Fig. 3(a) and (c-d), enlarging 
the collecting areas of the test detector can also enhance 
the quality of the images recovered by CS reconstruction. 
Fig. 2(b-d) and Fig. 3(b-c) also shows that CS recon- 
struction requires much shorter acquisition times than 
GI reconstruction for reconstructing ghost images with 
a good SNR, which have already been demonstrated in 
Ref. [3. 

In conclusion, ghost imaging via CS reconstruction is 



the result combining quantum optics with information 
theory. Based on CS reconstruction, except for allow- 
ing for shorter acquisition times for reconstructing ghost 
images with a good SNR, super-resolution ghost images 
can also be obtained even if the information from the 
object is partially collected and the best resolution of re- 
constructed images is closely related to the CCD camera 
registering the light field in the reference path. In practi- 
cal applications, enlarging the collecting areas of the test 
detector is much easier than making a lens with large ef- 
fective aperture. So ghost imaging via CS reconstruction 
is very useful to the imaging in the long wavelength ra- 
diation band and in the far field, microscopy, astronomy 
and so on. 
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